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Abstract 

By extending our recently proposed magnon-density-waves to low dimensions, we investigate, 
using a microscopic many-body approach, the longitudinal excitations of the quasi-one-dimensional 
(quasi-ld) and quasi-2d Heisenberg antiferromagnetic systems on a bipartite lattice with a general 
spin quantum number. We obtain the full energy spectrum of the longitudinal mode as a function 
of the coupling constants in the original lattice Hamiltonian and find that it always has a non- 
zero energy gap if the ground state has a long-range order and becomes gapless for the pure 
isotropic Id model. The numerical value of the minimum gap in our approximation agrees with 
that of a longitudinal mode observed in the quasi-ld antiferromagnetic compound KCUF3 at low 
temperature. It will be interesting to compare values of the energy spectrum at other momenta if 
their experimental results are available. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Ds, 75.50.Ee 
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I. INTRODUCTION 



The low temperature properties of many two-dimensional (2d) and three-dimensional 
(3d) quantum antiferromagnetic systems can be understood by Anderson's spin wave theory 
(SWT) and its extensions [l|, which provides correct description of the quantum corrections 
to the classical Neel states of the systems. For many purposes, the dynamics of these systems 
at low temperature can be considered as that of a dilute gas of weakly interacting spin-wave 
quasiparticles (magnons) with its density corresponding to the quantum correction to the 
classical Neel order; also present in these systems are the longitudinal fluctuations consisting 
of the multi-magnon continuum j^j]. 

On the other hand, due to the strong quantum fluctuations, the isotropic Id antiferro- 
magnets with low quantum spin numbers exhibit different low temperature properties, such 
as no Neel-like long-range order in the ground state and the quite different low-lying exci- 
tation states from the 2d and 3d counterparts. According to the exact solutions by Bethe 
ansatz, the natural low-lying excitation states of the Id spin-1/2 Heisenberg model have 
been shown corresponding to the spin-1/2 objects (spinons) which always appear in pairs, 
and the spin- wave-like excited states are actually the triplet states of a spinon pair |3J, con- 
trast to the doublet states from SWT. For the spin-1 Heisenberg chain, however, the triplet 
excitation states have a nonzero energy gap, first predicted by Haldane [4|. These theoret- 
ical predictions have been confirmed by the experimental results in the antiferromagnetic 
compound KCUF3 for the spin-1/2 chains [5'] and CsNiOs for the spin-1 chains 6]. 

Strictly, all real systems are three dimensional when temperature is low enough. The 
antiferromagnetic compounds KCuF 3 or CsNiCls are actually quasi-ld systems with very 
weak inter-chain couplings. In particular, the spin-spin couplings are ferromagnetic in the 
tetragonal basal planes of KCuF 3 and antiferromagnetic in the hexagonal planes of CSMCI3. 
Many parent compounds of the high-T c superconducting cuprates or the ion-based pnictides 
are also quasi-2d antiferromagnetic systems with very weak inter-plane couplings {7, 8]. 
Therefore, there is a 3d magnetic long-range order with a nonzero Neel temperature XV for 
all these systems and one expects SWT should provide a qualitatively correct description 
for some if not all the low-temperature dynamics of these quasi-ld or quasi-2d systems. One 
interesting question is whether or not some ld-type excitations, such as the longitudinal 
part of the triplet spin-wave excitations of the pure Id systems can survive in the ordered 
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phase when the inter-chain couplings are present. In fact, there is now ample evidence of 
the longitudinal excitation states in various quasi-ld structures with the Neel-like long-range 
order at low temperature, including the hexagonal ABX^type antiferromagnets with both 



spin quantum number s = 1 (CsNiCl 3 and RbNiCl 3 )_|9|, ll(| and s = 5/2 (CsMnI 3 ) \l2\ 
and the tetragonal structure of KCuF 3 with s = 1/2 13J. More recently, a longitudinal mode 
was also observed in the dimerized antiferromagnetic compound TlCuCl 3 under pressure with 



a long-range Neel order 



14|. To our knowledge, no observation of any longitudinal mode 



in the quasi-2d antiferromagnets has been reported yet. Clearly, such longitudinal modes, 
which correspond to the oscillations in the magnitude of the magnetic order parameter, are 
beyond the usual SWT which only predicts the transverse spin- wave excitations (magnons). 
There have been several theoretical investigations in these longitudinal modes, all using the 
field theory approach, sueh as a simplified version of Haldarre^ theory for the spin-1 systems 
[15| or the sine-Gordon theory for the spin-1/2 systems 16|, [17|, and both treating the inter- 
chain couplings as perturbation. A phenomenological field theory approach focusing on the 
spin frustrations of the hexagonal lattice of the ABX^-type antiferromagnetic systems has 
also been proposed [181 ] . Common to all these field theory approaches is the need to take 
the continuum limit with a number of fitting parameters. By proper choice for the values 
of the fitting parameters, general agreements with the experiments mentioned earlier have 
been found, although there are still some disagreements particularly for the data away from 
the minimum energy gap at the antiferromagnetic wavevector [l^ . 

We recently proposed a microscopic many-body theory for the longitudinal excitations of 
spin-s quantum antiferromagnetic systems, using the original spin lattice Hamiltonians 19]. 
The basic physics in our analysis is simple: by analogue to Feynmann's theory on the low- 



lying excited states of the helium-4 superfluid [20j, we identify the longitudinal excitation 
states in a quantum antiferromagnet with a Neel-like order as the collective modes of the 
magnon-density waves, which represent the fluctuations in the long-range order and are 
supported by the interactions between magnons. These longitudinal excitation states are 
constructed by the s z spin operators, contrast to the transverse spin operators s ± of the 
magnon states in Anderson's SWT. These modes are referred to as the collective modes 
of the magnon-density waves because of the fact that s z is the magnon-density operator 
in these systems. The energy spectra of these collective modes can be easily derived by a 
formula first employed by Feynmann for the famous phonon-roton spectrum of the helium 
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superfluid involving the structure factor of its ground state. Nevertheless, we now realize 
that the precise form for the definition of the longitudinal state in our earlier work is not 
quite correct and we have now slightly modified the definition and, indeed, we find the 
corresponding values of the energy spectra in an approximation using the SWT ground state 
are in general much lower than before. We extend our analysis to the Id model and find that 
in the isotropic limit the longitudinal mode has a gapless spectrum. Interestingly, this gapless 
spectrum in our approximation is degenerate with the doublet spin-wave spectrum of SWT, 
hence making it triplet, in good agreement with the exact triplet spin- wave spectrum of the 
spin-1/2 Heisenberg model [3]. The application of our analysis to the quasi-ld and quasi- 
2d systems is straightforward and hence more detailed comparison with the experiments 
is now possible. Our numerical results for the spin-1/2 quasi-ld compound KCuF 3 show 



the minimum gap value of the 
obtained from the experiments 



ongitudinal energy spectrum in agreement with the value 



131 ] . This is particularly satisfactory since our analysis has 



no fitting parameters except the coupling constants in the original lattice Hamiltonian. As 
our microscopic approach is able to obtain the full spectrum of the longitudinal mode, it 
will be interesting to compare the values at other regions if their experimental results are 
available. 

We present our general theory of the magnon-density-waves briefly in Sec. II, with nu- 
merical results calculated in an approximation using the SWT ground state for the simple 
cubic and square lattices and its extension to the Id models in Sec. III. We then discuss its 
application to quasi-ld and quasi-2d systems in Sec. IV, including the quasi-ld compound 
KCUF3. We summarize our work and discuss possible observations of the longitudinal modes 
in some quasi-2d systems in Sec. V. We also discuss the approximations employed in our 
analysis and their possible improvements in the last section. 



II. MAGNON-DENSITY-WAVES IN ANTIFERROMAGNETS 

We consider an antiferromagnetic system as described by an V-spin Hamiltonian H on a 
bipartite lattice. The classical ground state is given by the Neel state with two alternating 
sublattices a and b, where we assume the spins on the a-sublattice all point up in the in- 
direction of the spin space and the spins on the 6-sublattice all point down. This Neel state 
describes the perfect two-sublattice long-range order. In this article, we shall exclusively use 
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index i for the sites of the a-sublattice, index j for the sites of the 6-sublattice, and index / 
for both sublattices. The excited states are given by the spin-flipped states with respective 
to the Neel state and are commonly referred to as magnons, the quasiparticles of magnetic 
systems in general. 

The quantum ground state of H in general differs from the classical Neel state by 
a correction, the long-range order is hence reduced. For many purposes, as described by 
the spin- wave theory (SWT) this quantum correction in most 2d and 3d models is well 
represented by a gas of magnons whose density p directly gives the correction as 

(sZ) g = s-p, (1) 

where s is the spin quantum number, s z is the z-component of spin operator on the lattice site 
i, and the expectation (. . ) g is with respect to the ground state \^f g ). Similarly, (sf) fl = — s+p 
for the 6-sublattice with the same density p. Therefore, operators s z corresponds to the 
magnon-density operators, contrast to the spin-flip operators s ± which correspond to the 
magnon creation/destruction operators. Clearly, there are two types of the magnons due to 
the two sublattice structures. Anderson's SWT can be most simply formulated by bosonizing 
the two sets of these three spin operators, s z and s ± , on the two sublattices respectively. 
For example, the quantum correction to the classical Neel state by the linear SWT gives the 
magnon density of p = 0.078 per lattice site for the spin-1/2 Heisenberg model on a simple 
cubic lattice, and of p = 0.197 per lattice site for the same model but on a square lattice. 

Due to the interactions between the magnons, it may be necessary to consider the states 
of the magnon-density waves (MDW). These states may not be well defined in the 3d systems 
where the magnon density is very dilute and the long-range order is near perfect with little 
fluctuations. In the low dimensional systems, however, the magnon density may be high 
enough to support these longitudinal waves. In terms of microscopic many-body theory, 
these MDW states are the longitudinal excitation states constructed by applying the density 
operator s z on the ground state in a form as s z \^ g ), similar to Feynmann's theory of the 
phonon-roton excitation state of the helium superfluid, where the density operator is the 



usual particle density operator [20]. These longitudinal states may be compared to the 



quasiparticle magnon states which are constructed by the transverse spin-flip operators as 



s \^f g ). The above discussion underlines the main idea in our earlier papers 



19(| , whose main 



purpose is to outline a general framework for the excitation states of both quasiparticles 
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and quasiparticle-density waves for a general quantum many-body system in our variational 



coupled-cluster method 
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In more details, following Feynmann, the MDW excitation state with momentum q in an 
antiferromagnet is given by 

\* a q )=X a q \* g ), (2) 

where excitation operator X®, in the linear approximation, is the sublattice Fourier trans- 
formation of the magnon density operators si of the a-sublattice, 

x : = JlT,^- r '^, ?>o, (3) 

with the condition q > required because of its orthogonality to the ground state \^f g ) in 
which sf otal = 0. The excitation energy spectrum in this linear approximation can be derived 

as, 

N a (a) 

*•<*) - s# <«) 

where N a (q) is given by a double commutator, 



N a (q) = \([X%, [H, X$]) B , (5) 



and the state normalization integral S a (q) is in fact the structure factor of the a-sublattice, 

S a (q) = (X%X«) g = |E eiq ' (ri_ri,) «4) 9 - ( 6 ) 
Similarly, we have the MDW excitation state X b \^/ g ) with operator 

K = JlH^s z v g>0, (7) 

and the corresponding energy spectrum E b (q) for the 6-sublattice. Due to the lattice sym- 
metry, the spectra E g and E b q are degenerate. However, these two excitation states are not 
orthogonal to each other because of the couplings between the spins on the a-sublattice and 
the spins on the o-sublattice. We therefore need to consider their linear combinations as, 

= x t\* 9 ) = ^(*-± *;)i*,>, (8) 

for the coupled MDW states. The corresponding energy spectra is similarly given by 
E ± {q) = N ± (q)/ S ± (q) with N ± (q) and 5 ,± (g) given by similar equations as Eqs. (5) and 



(6) respectively using excitation operators instead of X®. It seems that we have two 
longitudinal modes. But these two states with the energy spectra E ± (q) are actually 
the same state, with one equal to another by a substitution for the wavevector q — y q + Q 
where Q is the antiferromagnetic wavevector of the system (e.g, Q = (n, it) for the 2d square 
lattice model). This can be easily seen as the excitation operator X^ are in fact nothing 
but the Fourier transformations of the (staggered) magnon density operators sf or (— l) l sf 
respectively. We therefore only need to consider one of them. We choose = \^ q ) with 
its energy spectrum E~(q) = E(q), and write 

where N(q) and S(q) are calculated by using X~ = X q of Eq. (8). We notice the slight 
difference between this definition of the MDW states of Eq. (8) and that in our earlier paper 

n 

|l9| where we used the total density operator as hi = h(2s — sf + ~J2n s i+ n ) with z as 
the coordination number and n as the nearest-neighbor index. We now realize the use of 
operator hi (or its equivalent form, sf — - J2 n s t+n) ls n °t quite correct. Our current definition 
of Eq. (8) seems more natural as discussed in details above. Indeed, as we will see later, the 
values of the energy spectrum of the states defined by Eq. (8) in our approximation are in 
general much lower than before, with the maximum energy values about half of those of the 



earlier results [19 |. 



So far, in the above general analysis for the longitudinal excited states, the exact ground 
state \ty g ) is used for the ground state expectation values. The only approximation comes 
from the choice of the linear form in the excitation operators of Eqs. (3) and (7), and is often 
referred to as the single-mode approximation as viewed from the general expression of the 
dynamic structure factor. In the case of the helium superfluid, the double commutator can 
be simply evaluated as N(q) oc q 2 , and Feynmann [3] used the experimental results for the 
structure factor with S(q) oc q as q — y and hence derived the low-lying phonon spectrum 
E(q) oc q and the gapped roton spectrum around the peak of S(q). Jackson and Feenberg, 
however, used the variational results calculated from the Jastrow-type wavefunctions and 



obtained similar results [2l{]. In our earlier papers |l9j . we have demonstrated that these 
equations remain valid when the exact ground state \^> g ) is replaced by a variational state 
|^o) ail d furthermore, in the case of the quantum antiferromagnets as discussed here, our 
variational ground state |^ ) by the so-called variational coupled-cluster method in a first 



order approximation reduces to that of Anderson's SWT [22J]. Therefore, to this first order 
approximation which is what we focus on here, we apply the SWT ground state |\1> SW ) in 
all of our following calculations. We like to emphasize that SWT itself in its usual form 
cannot produce the longitudinal MDW excitations discussed here. We will discuss this 
approximation and its possible improvement in the last section. We present our numerical 
results of the MDW spectra E(q) for several models in the following two sections: Sec. Ill 
contains the results for the antiferromagnetic models on the simple lattices, while Sec. IV 
contains results for the more physical quasi- Id and quasi-2d systems. 

III. RESULTS OF MAGNON-DENSITY-WAVE SPECTRA IN SIMPLE LAT- 
TICES 

A. The spin-s XXZ Heisenberg model 

In this section, we present the numerical results for the energy spectra of the MDW states 
as discussed in the earlier section for the spin-s XXZ Heisenberg model on a simple cubic 
lattice and a square lattice. We then present the results for the Id model and discuss the 
convergent results in its isotropic limit. 

The spin-s XXZ Heisenberg model on a bipartite lattice is given by 



where the coupling parameter J > 0, index i runs over all a-sublattice only, index n runs over 
the z nearest-neighbor sites, and A (> 1) is the anisotropy parameter. The usual isotropic 
Heisenberg model is given by A = 1. The purpose of introducing the anisotropy is twofold: 
it is interesting on its own right and it also provides a way to obtain convergent results for 
the Id case in the isotropic limit as we will see later. 

Using the usual spin commutation relations, it is straightforward to derive the following 
double commutator as, 



i,n 



+ S 




(10) 



N(q) = i([X_ g , [H, X g }}) 9 = - Z J-(l +lq )(sts- +l ) 9 , 



(11) 



where 7 9 is defined as usual, 



(12) 



n 
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with the coordination number z, and (sfs^ +1 ) g is independent of the index % due to the lattice 
translational symmetry. The general expression for the structure factor S(q) contains an 
additional cross term compared to the sublattice counterpart S a (q) as 



S(q) = (XVQ. = S a (q) + -J2e i «< Ti - r > ) (s*s*) g . (13) 

Before we discuss any approximation, we notice that the double commutator in general 
behaves as, near the antiferromagnetic wavevector Q, 

iV(|q+Q|)ocg 2 , g^O, (14) 



similar to that of the helium superffuid 20 ]. 



Now we need a specific approximation for the ground state in order to evaluate 
the spin correlation functions (sfsj) g , (sfsf,) g , and (sfsf) g . As mentioned earlier, in this 
article we use as our first-order approximation the spin- wave ground state, |\l/ sw ), for these 
calculations. After defining the transverse spin correlation function g(r), 

9r = ^(4 s 7+r)gi (15) 

we derive the following results for its Fourier transformation, 

1 iJA 



and the sublattice structure factor, 



(16) 



11 1 

with the magnon density p = J2 q Pq- And, finally, the full-lattice (staggered) structure factor 
is given by, 

S(q) = S a (q) + Y,9 q ,g q - q , (18) 

We notice that in deriving the expressions of Eqs. (17) and (18) for the structure factors, the 
values for q = are excluded due to the condition q > in the definition of X q from Eqs. (3) 
and (7). Furthermore, the integrals in the structure factor involving function r j q i r ) q - q i clearly 
indicate the couplings between magnons. In all these formulas, the summation over q is 
given by 



?=(2^/-/ 9 ' (19) 



where d is the dimensionality of the system. The energy spectrum E(q) of Eq. (9) is obtained 
by calculating the values for N(q) and S(q) from the approximations of Eqs. (16-18). This 
longitudinal spectrum E(q) can be compared with the following transverse spin- wave spectra 
of the linear SWT [l|, 

E sw (q) = szJA^Jl - rf/A 2 . (20) 
In the following subsections, we present numerical results using the above approximations. 



B. Results for the simple cubic and square lattices 

We first consider the isotropic case A = 1 for the simple cubic lattice model for which gx is 
calculated as gi = (sfs~ +1 ) g ~ —0.13. The numerical values for E(q) near q — » are similar 
as given earlier [191 ]. with a large energy gap of about 0.99szJ at q = 0. But at other values 
of q, the energies are much smaller than before due to the different definitions of the density 
operator of Eq. (8) |23j. At the antiferromagnetic wavevector (AFWV) Q = (tt,tt, 7r), the 
spectrum has a larger gap of lAOszJ. As discussed before, this high energy 3d longitudinal 
mode may not be well defined and distinguishable from the multimagnon continuum. 

For the square lattice model at the isotropic point A = 1, g\ « —0.28. Similar to the 
earlier results [3], E(q) becomes gapless at both AFWV Q = (it, it) and q — > 0, due to the 
logarithmic behaviors from the structure factors (e.g., S(q) — > — lng hence E(q) oc — 1/lng 
as q — > 0). However, as discussed earlier, this logarithmic gapless spectrum of the square 
lattice model in fact is quite "hard" in the sense that any finite-size effect, anisotropy or 
interplane coupling to be discussed later, however small, will make a nonzero gap. For 
example, we consider a tiny anisotropy here with a value A — 1 = 1.5 x 10 -4 , which in fact is 
a typical value for the parent compound of the high-T c superconducting cuprate, La2Cu04 



24j , we obtain in our approximation the gap values at E(Q) ~ 0.7QszJ and E(q) ~ OAAszJ 
as q — > 0, both much larger than the corresponding magnon gap value of 0.02szJ from 
Eq. (20). We plot part of the spectrum with this anisotropy in Fig. 1, together with the 
spin-wave spectrum for comparison. The energy values at the two particular momenta 
(7r/2,7r/2) and (vr, 0) deserve attention, where 7 g = and the spin-wave spectrum gives the 
same value of szJ. The longitudinal spectrum E(q) at these two point has slightly different 
values, 1.36szJ and lAOszJ respectively. The energy difference at these two momenta has 
been used to indicate nonlinear effects due to magnon-magnon interactions in the more 
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accurate calculations for the isotropic Heisenberg model [25]. It is interesting to note that 
our longitudinal mode also show this difference. 



C. Results for the Id model 

We next consider the Id case. The SWT results in general for the isotropic Id case 
are not reliable as most integrals suffer from the well-known infrared divergence, e.g., the 
magnon density p — > oo as A — > 1, an unphysical result. Nevertheless, the value of the 
spin-wave spectrum of Eq. (20) is not far off that of the exact result by Bethe ansatz for 
the spin-1/2 model despite the different degeneracies (i.e., the spin-wave spectrum is doublet 
while the exact spectrum is triplet). The infrared divergence of the spin- wave results also 
occurs for the parameter cji in the numerator of the energy spectra in Eq. (9). We examine 
the behaviors of each integral in N(q) and S(q) in the isotropic limit A — > 1 and find that 
they all have the similar infrared divergence. For example, by numerical calculations, we 
find that 

gt oc --^-la(A- 1), asi^l, (21) 



agrees with the analytical results using the elliptical formula 26]. Furthermore, in the limit 
q — > 0, both S a (q) and S(q) behave as 

1 \n(A — 1) 

S{q) ->• ) ' , as q ->• and A ->• 1. (22) 

27T y/ A — 1 

Since the divergences in the numerator N(q) and the denominators S(q) precisely cancel out, 
we obtain finite results for the energy spectrum E(q) for the isotropic Id model. Interestingly, 
we find that these numerical values of E(q) coincide precisely with those of the linear spin- 
wave spectra of Eq. (20) for all values of q in the isotropic limit A — > 1. Therefore, our 
longitudinal spectrum and the doublet transverse spin- wave spectrum constitute a triplet, 
in good agreement with the following exact triplet spectrum for the spin-1/2 model by Bethe 
ansatz first derived by des Cloizeaux and Pearson 3|], 

E( q ) = ^L sin q. (23) 

The different factor irJ/2 of the above exact result comparing to the value of J by the linear 
SWT of Eq. (20) with z = 2 clearly comes from the nonlinear effects beyond our simple 
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approximation employed here. We also notice that our analysis here in the approximations 
employed is not able to produce the Haldane gap for the isotropic spin-1 chain. 

For the anisotropic Id model (i.e., A > 1), the triplet spectra split and the values of the 
longitudinal spectrum E(q) are larger than those of the doublet spin-wave spectrum, similar 
to the cases of the 2d and 3d models discussed earlier. We plot this E(q) for A = 1.1 in 
Fig. 2 as an example. The gaps for E(q) are about l.lGszJ and l.QAszJ at q = and it 
respectively, comparing with OAQszJ of the spin-wave spectrum at both points. 



IV. MAGNON-DENSITY WAVES IN QUASI-1D AND QUASI-2D SYSTEMS 

A. Quasi-ld and quasi-2d antiferromagnets on bipartite lattices 

A generic quasi-ld and quasi-2d antiferromagnetic Hamiltonian on a bipartite lattice is 
given by, 

H = J Sj • s i+ni + Jj_ Sj ■ s i+n2 , (24) 

where index i as before runs over all a-sublattice sites with index ri\ over the nearest- 
neighbor sites along the chains and HI2 over the nearest-neighbor sites on the basal planes, 
J is the coupling constant along the chains and J± is the counterpart on the basal planes. 
We consider the model with both J and J± > 0. The quasi-ld model corresponds to the 
case of J±/ J <C 1, the quasi-2d model to the case of J±/J ^> 1, and the 3d model is given 
by J± = J. This Hamiltonian has been studied for the case of the quasi-ld systems with 
J±/J <C 1 by SWT 0, Q- I n particular, the SWT ground state was used to evaluate the 
corrections due to the kinematic interactions to the order parameter p. The longitudinal 
modes were not discussed. 

All of our earlier formulas for the longitudinal mode at the isotropic point A = 1 remain 
the same after the following replacements 

2 

z -> z' = 2(1 + 2£), 7 9 ->■ 7^ = — [cosq z + £(cosq x + cosq y )} , (25) 

where £ = Jj_/J. This is true also for the spin-wave spectrum E sw (q) of Eq. (20). We 
notice that the spin-wave spectrum is gapless at zone boundaries, the longitudinal mode 
E(q) of Eq. (9) however has nonzero gaps for any £ > 0, at which there is a long-range 



order 129, 



27| . In Fig. 3, we present our results for the spectrum, denoted as E qld (q), of 
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the quasi-ld model with A = 1 and £ = 1.05 as an example, together with the spin- wave 
spectrum E sw (q). The gaps for E qld (q) at q — > and Q = (7r,7r,7r) are 0.78sz'J and 
1.21sz'J respectively. Fig. 3 also includes our results for a quasi-2d model with A = 1 and 
l/£ = J/Jj_ = 10~ 3 , denoted as i? q2d (g). The gap values for this quasi-2d spectrum at g — > 
and Q = (tt, it, 7r) are about 0.47s;?' J and 0.80S2;' J respectively. We also notice that at the 
particular two momenta (vr/2, tt/2, 0) and (it, 0, 0), where the linear spin-wave spectrum has 
the same value of sz'J but the longitudinal mode has slightly different values, 1.36sz'J and 
lAlsz'J respectively, due to magnon-magnon interactions as discussed earlier. This quasi-2d 
model may be relevant to the parent compounds of the high-T c cuprates, where the effective 
interlayer couplings J/Jj_ between the Cu0 2 planes are estimated to be between 10~ 2 and 
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28|. 



B. Quasi-ld model with KC11F3 structure 

In the experimentally well studied quasi-ld compound KCuF 3 , the strong spin couplings 
along the chains are antiferromagnetic but the weak couplings on the basal plane are ferro- 
magnetic. This compound can be described by the following Hamiltonian model, 

H = J Y, S L- Sia+m -Y\^ Sla ' Sla+n2 + ^ Sl *' s ^+" 2 ) ' ( 26 ) 

la,ni \la,ri2 h,ri2 / 

whose classical Neel state consists of two alternating planes, with all the spins on the a- 
plane pointing up and labeled by index l a and all the spins on the 6-plane pointing down 
and labeled by index l^. In Eq. (26), the nearest-neighbor indices ni and n 2 are as defined 
before with ri\ along the chains and n 2 on the basal planes, and both J and J± > 0. The 
spin-wave spectrum is derived as 

E sw (q) = 2sJ^/Al-cos 2 q z , (27) 

where A q is defined as 

A 9 = l + 2£(l- 7 f), e = y, (28) 

with 7 2d = (cos q x + cos q y ) /2. It is easy to check that when J = 0, we recover the spin-wave 
spectrum of the 2d ferromagnetic model and that when J± = 0, we recover the spin-wave 
spectrum of the Id antiferromagnetic model. For the longitudinal energy spectrum E(q) of 
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Eq. (9), the double commutator is now given by a different form as 

N(q) = -J(l + cos q z )(stsT a+ni ) g + 2J ± (1 - 7 f )«^ +n2 )„ (29) 

and the structure factor S'(g) is as given before by Eqs. (6) and (13) in general forms and 
by Eqs. (17) and (18) in our approximation using the similar SWT ground state with the 
anisotropy parameter A = 1. We notice that in Eq. (29), the two spin operators in the first 
correlation function (sf a sj' +ni ) g are from the two sublattices respectively as before, but in the 
second correlation function ( s f a sT a +n 2 )g> they are from the same sublattice. So we still name 
the first one as before by g\ = (s^"s^ + ) g /2s but the second one as g[ = (si a sj a+n2 ) g /2s. 
Using the similar SWT ground state, we obtain, for their Fourier transformations, 

1 cosq z 

9 « = ~9 m 2 (30) 

Z ,/A„- cos z q z 



i 

and 



9 * = p <> = 2 /A2 2 ( > 



respectively. We notice the quite different expressions for g q and g' as expected. We present 
our numerical results for E(q) in Fig. 4, together with E sw (q) of Eq. (27) for comparison, 
using the experimental values for the coupling constants, J ~ 34 meV, J± ~ 1.6 meV and 
s = 1/2 [13]. Different to the longitudinal modes in other systems discussed earlier, we 
find that E(q) has a smaller gap of about 0.63J^ 21.4 meV at AFWV Q = (0,0, tt), and 
a larger gap of about 0.85J~ 28.9 meV at q — > 0. This gap value of 21.4 meV at AFWV 
is about 40% higher than the experimental value of about 15 meV. The field theory by 
Essler et al. produces a gap value of about 17.4 meV 17j. However, there is uncertainty 



in the estimate value of the interchain coupling constant J±. Lake et al. seem to have 
used the theoretical formula Eq. (56) in Ref. [lTj to obtain J± = 1.6 meV = 0.047 J. By 



different methods 



aollaii 



J± was estimated to be 0.01J ~ 0.016J. Using this estimate of 



£ = J±/J = 0.01, we obtain the minimum gap value of 11.9 meV at AFWV and 16.8 meV 



at q — 0. Naively, if we choose about the midpoint between the values of Refs. [13j and 30J, 
J± = 0.85 meV with £ rj 0.025, we obtain the minimum gap value of 0.49J = 16.8 meV 
at AFWV and 0.68J = 23.2 meV at q — 0, in good agreement with the experiment for the 
minimum gap [l3|. Furthermore, with this value of J± = 0.85 meV, the linear spin-wave 
spectrum gap at q = (7r,0,7r) is E sw (q) rs 0.32J = 10.9 meV, very close to the gap value 
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of 11 ± 0.5 meV by the experiment [30j. The longitudinal mode E(q) is nearly flat in the 
region (r],0, 7r) ~ (77, 0, 7r) with rr < rj < 0, with the gap value about 0.59J = 20.1 meV at 
(7r, 0, 7r). It will be very interesting indeed to compare with experimental results if available 
for the whole spectrum. 



V. SUMMARY AND DISCUSSION 



In summary, we have investigated the longitudinal excitations of various quantum anti- 
ferromagnets based on our recently proposed magnon-density-waves. Our numerical results 
show that the longitudinal mode always has a nonzero gap so long the system has a Neel- 
type long-range order and becomes gapless in in the limit of the Id isotropic model. In 
particular, the spectrum of the longitudinal mode in our approximation is degenerate with 
the doublet spin-wave spectrum of SWT in the limit of the isotropic Id model, in agreement 
with the triplet spin- wave spectrum of exact results for the spin- 1/2 model by Bethe ansatz 
j^]. In the case of the simple cubic lattice model, the longitudinal mode with high energy 
values may not be well defined since there is little fluctuations in the nearly perfect classical 
long-range order. In the quasi-ld and quasi-2d models, where the quantum correction is 
large and the magnon density is significant, the magnon- density waves may be observable. 
Indeed, there are now ample evidence of the longitudinal modes in several quasi-ld systems 
as mentioned earlier in Sec. I. In particular, for the quasi-ld compound KCuF 3 , our value 
for the minimum gap is in agreement with the experimental value [ljj. It will be interest- 
ing if more experimental results for the spectrum away from the minimum are available for 
comparison. 

It is also interesting to note that the longitudinal modes were observed in the ABX 3 -type 
antiferromagnets with both s = 1 [9, [h3| and s = 5/2 [ll, 12], clearly indicating that the 
modes are more general in their physics, independent of the mechanism which generates 
Haldane gap of the Id model. The phenomenological field theory model with five fitting 
parameters employed by Affleck is derived from Haldane's theory of the spin-1 chain 151 ] . It 
will be interesting to apply our general microscopic analysis presented here to the ABX^-type 
antiferromagnets where the basal plane is hexagonal and the corresponding Neel-like state 
has three sublattices rather than two sublattices discussed here. Other systems where we 
can apply our analysis for the magnon-density-waves include the quasi-2d systems where the 
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next-nearest-neighbor antiferromagnetic couplings, in addition to the usual nearest-neighbor 
couplings, are present. These additional couplings cause quantum frustrations and the Neel- 
like order is further reduced hence greater the magnon density to support the magnon- density 
waves. Of particular current interest is the the parent compounds of the newly discovered 
high-T c superconducting ion-based pnictides where such next-nearest-couplings are believed 
to be significant [29 ]. 

Finally, we want to point out that there are two major approximations in our analysis 
here. The first is the linear operators X q employed in constructing the excitation states 
and the second is the SWT ground state |\I/ SW ) employed in evaluating all the correlation 
functions involved. In regard to the first approximation, it is interesting to consider the 
case of the phonon-roton spectrum of the helium superfluid, where after inclusion of the 
nonlinear terms due to the couplings to the low- lying phonons (i.e., the so-called backflow 
correction), the values of the roton gap are reduced by about half to near the experimental 



values 



20 



211- 



Clearly, the effects due to the couplings between the longitudinal modes and 
the gapless magnons in the antiferromagnetic systems also deserve further investigation. In 
regard to the second approximation, i.e., the SWT ground state employ in our calculations, 
improvement can be obtained by using better ground state functions available by more 
sophisticated microscopic many-body theories such as the coupled-cluster method [22|, [32 ] 
and, particularly, its most recent extension where the strong correlations are included by 
a Jastrow correlation factor 33j. We believe the quasi-ld and quasi-2d antiferromagnetic 
systems as studied here are good theoretical models from both the view point of the field 
theory approach which deal with most effectively the nonlinear effects of the Id systems 



1* and of the microscopic many-body theory approach wh.eh provides genera., systematic 
techniques in dealing with many-body correlations in plethora of quantum systems [34|]. The 
two theoretical approaches complement one another in study of these models and we wish 
to report our progress in these investigations in near future. 
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FIG. 1: The energy spectrum E(q) for the longitudinal mode of Eq. (9) for the square-lattice 
Heisenberg model of Eq. (10) with an anisotropy A — 1 = 1.5 X 10~ 4 , together with the linear 
spin- wave spectrum E sw (q) of Eq. (20) for comparison. This tiny anisotropy is a typical value for 
the parent compound of the high-T c superconducting cuprate, La2Cu04. 
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FIG. 2: Similar to Fig. 1 but for the Id model with the anisotropy A = 1.1. In the isotropic limit 
of A = 1, E{q) approaches to E sw (q), forming a triplet spectrum as described in details in text. 
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FIG. 3: Similar to the earlier figures but for the quasi-ld and 2d systems of Eq. (24), with parameter 
£ = J±/J = 0.01 for the quasi-ld spectrum E q - ld (q) and £ = 10 3 for the quasi- 2d spectrum E<& d (q). 
The spin- wave spectrum E sw (q) is for the quasi-ld model. 
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FIG. 4: Similar to Fig. 3 but for the quasi-ld structure of KCUF3 as described by Hamiltonian of 
Eq. (26), with parameter £ = J±/J = 1.6/34 from the experiment 
E sw (q) is given by Eq. (27). 



13l | . The spin- wave spectrum 
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